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Abstract 

Osteonecrosis of the jaw, in association with bisphosphonates (BRONJ) used for treating osteoporosis or cancer, is a severe 
and most often irreversible side effect whose underlying pathophysiological mechanisms remain largely unknown. 
Osteocytes are involved in bone remodeling and mineralization where they orchestrate the delicate equilibrium between 
osteoclast and osteoblast activity and through the active process called osteocytic osteolysis. Here, we hypothesized that (i) 
changes of the mineralized tissue matrix play a substantial role in the pathogenesis of BRONJ, and (ii) the osteocyte lacunar 
morphology is altered in BRONJ. Synchrotron |xCT with phase contrast is an appropriate tool for assessing both the 3D 
morphology of the osteocyte lacunae and the bone matrix mass density. Here, we used this technique to investigate the 
mass density distribution and 3D osteocyte lacunar properties at the sub-micrometer scale in human bone samples from 
the jaw, femur and tibia. First, we compared healthy human jaw bone to human tibia and femur in order to assess the 
specific differences and address potential explanations of why the jaw bone is exclusively targeted by the necrosis as a side 
effect of BP treatment. Second, we investigated the differences between BRONJ and control jaw bone samples to detect 
potential differences which could aid an improved understanding of the course of BRONJ. We found that the apparent mass 
density of jaw bone was significantly smaller compared to that of tibia, consistent with a higher bone turnover in the jaw 
bone. The variance of the lacunar volume distribution was significantly different depending on the anatomical site. The 
comparison between BRONJ and control jaw specimens revealed no significant increase in mineralization after BP. We 
found a significant decrease in osteocyte-lacunar density in the BRONJ group compared to the control jaw. Interestingly, the 
osteocyte-lacunar volume distribution was not altered after BP treatment. 
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Introduction 

Under healthy conditions bone undergoes continuous remod- 
eling to adapt to spatially and temporarily variable demands, 
through a delicate equUibritim between resorption and formation, 
which is performed by osteoclast and osteoblast cells, respectively. 

Bisphosphonates (BP), which are commonly prescribed in the 
treatment of osteoporosis and bone metastasis, have been shown to 
reduce significantly the risk of fracture [1,2]. The action of BP 
relies on the reduction of bone resorption by inhibiting osteoclast 
activity. However, a severe and most often irreversible adverse 
effect of high-dosage BP treatment is the potential occurrence of 
osteonecrosis of the jaw [3-7]. Although multiple hypotheses have 
been formulated recently, the underlying pathophysiological 



mechanisms of bisphosphonate-related osteonecrosis of the jaw 
(BRONJ) are stiU not completely understood [8-1 1]. 

Bone remodeling results in a heterogeneous distribution of 
mineralized tissue units, with variable degrees of mineralization. 
This heterogeneity can be assessed from the bone mineralization 
density distribution by using techniques such as quantitative 
backscattered electron imaging (qBEI) [12-14], microradiography 
[15], or synchrotron radiation micro-CT [14,16]. The degree of 
mineralization of bone is a quality factor that influences tiie 
mechanical properties of bone [17]. 

Osteoclast and osteoblast activity is thought to be orchestrated 
by osteocytes, which are the most abundant type of bone cell and 
form a well-distributed network within the mineralized matrix 
[18]. These cells reside in cavities called lacunae measuring several 
hundreds of |J.m' in volume, and are interconnected through cell 
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Table 1. Sample 


details such 


as age, gender, anatomical reg 


on, BP treatment and duration of the BP treatment are listed. 




Intcrndl ssmple 










Duration of BP 




name 


Site 


Gender 


Region 


Age 


BP treatment treatment 


Underlying diagnosis 


tib29 


tibia 


male 


midshaft 


29 






tib56 


tibia 


male 


midshaft 


56 






tib88 


tibia 


male 


midshaft 


88 






fem5RF68 


femur 


female 


midshaft 


68 






feml5RF66 


femur 


female 


midshaft 


66 






fem11LF64 


femur 


female 


midshaft 


64 






femllLF87 


femur 


female 


midshaft 


87 






fem1LF70 


femur 


female 


midshaft 


70 






fem2RM60 


femur 


male 


midshaft 


60 






fem1LM71 


femur 


male 


midshaft 


71 






jaw3wk5 


jaw (BRONJ) 


female 


37 


n. k. 


Z >1 year 


n. k. 


jaw! wk4 


jaw (BRONJ) 


female 


35 


75 


Z 13 months 


Mammary-carcinoma 


jawlwkSB 


jaw (BRONJ) 


female 


n. k. 


70 


Z 16 months 


Plasmacytom 


jaw2wk2A 


jaw (BRONJ) 


female 


45 


70 


A 17 months 


Osteoporosis 


jawBwklA 


jaw (BRONJ) 


female 


13 


72 


A 1 0 years 


Osteoporosis 


jaw2wk36A 


jaw (BRONJ) 


female 


36/37 


74 


A 19 months 


Osteoporosis 


jaw2wk1 


jaw (BRONJ) 


female 


15 


44 


Z 1 year 


Mammary-Carcinoma 


jaw3wk4 


jaw (BRONJ) 


female 


45 


n. k. 


Z 2 years 


n. k. 


jaw1mk4 


jaw (BRONJ) 


male 


15 


84 


Z 1 0 years 


Prostate-carcinoma 


jaw2mk1 


jaw (BRONJ) 


male 


16/17 


81 


Z 14 months 


Multiple myeloma 


jaw2mg2 


jaw (control) 


male 


36 


44 






jaw2mg6A 


jaw (control) 


male 


48 


19 






jaw2mg3 


jaw (control) 


male 


48 


54 


- 




jaw2mg4 


jaw (control) 


male 


47 


42 






jaw3mg6 


jaw (control) 


male 


33/34 


27 






jaw2wg5 


jaw (control) 


female 


42 


42 






jaw2wg3B 


jaw (control) 


female 


17 


40 






jaw2wg1 


jaw (control) 


female 


16 


68 






jaw2wg4 


jaw (control) 


female 


42 


47 






'Z': intravenous administration of Zoledronate, 'A': oral administration of Alendronate, 


'n.k/: not known. 
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dendrites extending in tiny canals called canaliculi. The canalic- 
ular diameter of human bone has been reported to be within the 
range between 200 and 900 nm [19]. 

Osteocyte activity is thought to be stimulated by biological and 
mechanical signals [20] . The morphology of the lacuno-canalic- 
ular network (LCN) is believed to be related to the mechan- 
osensation and mechanotransduction processes of osteocytes [21- 
3 1] . Furthermore, the LCN ensures the transport of cellular waste 
and nutrients [23] . The LCN has also been reported to be essential 
for micro-crack repair by triggering bone remodeling [32]. In 
addition to their mechanical function, it is hypothesized that 
osteocytes regulate mineral metabolism, e.g. bone phosphate 
metabolism [33,34]. 

It was recently shown in a murine lactating model that not only 
osteoclasts are able to resorb bone matrix, but also that osteocytes 
remodel their peri-lacunar and peri-canalicular matrix [35]. 
Alterations in lacunar size have also been observed in response 
to changes of the mechanical environment, for example enlarged 
lacunae were reported in mice following space flight [36], or after 
glucocorticoid treatment [37]. In ovariectomized rats, both 



lacunar size and density were found to be altered in newly-formed 
bone after antiresorptive and anabolic pharmaceutical treatment 
[38]. 

So far, investigation of the three-dimensional (3D) structure of 
the LCN has been limited by the imaging techniques available 
[21,39]. Synchrotron radiation micro-computed tomography (SR 
|tCT) enables 3D imaging of bone tissue at the cellular length scale 
and has been shown to be an appropriate tool for investigating 3D 
lacunar morphology [40-43]. At the sub-micron resolution, SR 
PlGT enables 3D imaging of the LCN with a large field of view 
[44,45]. Synchrotron X-ray nano-CT with phase contrast, which 
provides sensitivity to the mass density variations that is several 
orders of magnitude higher than conventional attenuation contrast 
SR|j.CT, was recently used to investigate the bone LCN and the 
3D coUagen orientation at the nanometer length-scale [46,47] . 

However, very limited data is currently available at the sub- 
micron length scale for human jaw bone, both in terms of the 
distribution of osteocyte lacunae and in terms of mass density 
distribution [48]. 
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A B 




Figure 1. Slides of a reconstructed volume corresponding to a control jaw specimen are shown in A (x-y-plane) and B (x-z-plane). 

The white dashed lines in A and B indicate where A is located in B and vice versa. C shows a minimum intensity projection, the projection range is 30 
pixels (10.5 |j,m). D shows the lacunae mask corresponding to C, also in the form of an intensity projection (z-range = 30 pixels). Color bar in mass 
density of a.u. 

doi:1 0.1 371 /journal.pone.0088481 .gOOl 
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Figure 2. In (A) a relative mass density distribution is shown. The peak at around -100 is due to the air outside the sample, and the broad 
peak at about 200 is due to the glue used to attach the sample that managed to travel up the sample. In (B) the MDD"^ after segmentation is shown, 
and the parameters derived from it are indicated. 
doi:1 0.1 371/journal.pone.0088481 .g002 



PLOS ONE I www.plosone.org 



3 



February 2014 | Volume 9 | Issue 2 | e88481 



Mass Density and Lacunar Properties in Bronj 



In the present study, we investigated the difTerences in lacunar 
morphology and peri-lacunar tissue properties at the sub- 
micrometer length scale in human jaw bone tissue samples 
obtained from both healthy subjects and patients suffering fi'om 
BRONJ, Additionally, cortical bone samples collected from the 
femur and tibia of donors not treated with BP were analyzed and 
compared to BRONJ and healthy control jaw bone samples. AU 
investigations were based on SR [iGT using phase contrast images 
measured with a 350 nm isotropic pixel size. Coupled with phase 
retrieval, this simultaneously provides information about the 3D 
distribution of the osteocyte lacunae and the local mass density of 
bone [46]. In order to quantify the mass density distribution, 
which is closely related to the bone mineral density distribution 
(BMDD), we used indices similar to those introduced for BMDD 
[13]. We hypothesized that (i) the bone turnover of the healthy 
human jaw bone, assessed by analyzing the mass density 
distribution, would be increased in comparison to the other 
anatomical sites and (ii) that extracellular matrix density and 
lacunar volumes of samples originating from patients suffering 
from BRONJ would be altered in order to compensate for the 
mineral homeostasis disturbed by the inhibited osteoclast activity. 
We therefore reported and compared the osteocyte lacunae 
volume distribution and the spatial arrangement of lacunae, as 
well as descriptors for the mass density distribution of the peri- 
lacunar tissue. 

Materials and Methods 

Ethics Statement 

Ethical approval for the jaw samples was granted by the 
Arztekammer Bremen (Studien-Nr. 310). AU donors signed an 
informed consent form. The present study also involved the use of 
cadaver specimens. 

Ethical approval for the femur samples was granted by the 
Ethical Commission of the Medi[:al University of Vienna, see [49]. 
Ethical approval for the tibia was granted by the University of 
California, San Francisco Committee on Human Research, see 
[50]. 

Specimen Preparation 

Nineteen human jaw bone sections (blocks of about 1-3 mm'^ in 
size) were extracted from 12 female and 7 male donors, of whom 8 
female and 2 male donors were suffering from BRONJ. The 
healthy control samples were from debris obtained during tooth 
removal. The BRONJ samples were obtained from surgeries 
necessary for the treatment of the necrosis. Furthermore, 7 
cadaver specimens originating from the human femoral midshaft 
and 3 cadaver specimens originating from the human tibia 
midshaft of other donors were included in the present study. 
Detailed information on g(^ndcr, donor age, anatomical origin, 
and BP treatment can be found in I'aljle 1 . 

The femora were collected and prepared as described previously 
[49]. After removal, the jaw bone sections were embedded in 
Tissue-Tek® O.C.T.™ (Sakura Fintec Europe B.V., Alph en aan 
den Rijn, Netherlands) solution and stored frozen at — 20°C until 
further processing. Following thawing, the specimens were drilled 
down to a diameter of about 0.5 mm using a high precision lathe 
[51]. For the BRONJ samples, this sub-volume was selected from a 
region in which no necrotic tissue had been observed. The cut 
bone samples were fixed in 70% ethanol for transport. Tibia bone 
cores were machined with a coring tool and precision circular saw 
to a length of 4 mm and diameter of 4 mm, similar to the sample 
preparation described in [50] . Then, to fit the field of view of the 
imaging setup, the tibia samples were cut from cortical bone 



(diameter: 500 nm, height: 1 mm) using a high precision drilling 
machine. About 12 hours before imaging, the samples were placed 
inside the measurement hutch in order to allow adjustment to 
humidity and temperature. 

Synchrotron Radiation Phase Contrast |j,CT 

The SR ^.CT data were obtained at ESRF (European 
Synchrotron Radiation Facility, Grenoble, France) at beamline 
ID22NL The X-ray beam was focused using Kirkpatrick-Baez 
reflective optics [46]. The scans used in the present study were 
performed for each sample by collecting 1201 projections, each 
with 0.2 s exposure time, over a total range of 360°. The energy 
was set to 16.874 keV, and the sample-detector distance was 
282 mm, resulting in a (350 nm)'' isotropic voxel size in the 
reconstructed image. Due to the coherence of the synchrotron 
source, the intensity of the recorded radiograph includes phase 
contrast [52,53]. Reconstruction was performed using Paganin's 
method [54], coupled to the conventional filtered back projection 
algorithm. 

In the Paganin method, the phase is retrieved by simply 
assuming a linear relationship between the absorption index (P) 
and the refractive index decrement (5). For cortical bone (ICRU- 
44), the S/p (delta/beta) ratio at the given energy was set to 199 
based on the XOP software [55]. The high ratio of delta/beta 
demonstrates the higher sensitivity for imaging the phase (8) 
compared to imaging the attenuation (P). The reconstructed 3D 
image made of 2048'' voxels corresponds to a map of the refractive 
indices stored in units of 2Jl/X,, with X being the wavelength of the 
X-ray beam (here X, = 0.0735 nm). This map is linearly related to 
mass density [46] which was shown to be associated with the 
degree of mineralization [56]. The spatial resolution with these 
settings allowed an easy distinction of osteocyte lacunae and larger 
pores from the mineralized tissue matrix (Fig. lA-C). However, 
the canalicular network could be not resolved. 

Image Segmentation 

In order to segment osteocyte lacunae inside the bone tissue 
volume (BV), the histogram of the whole 3D image was computed 
(Fig. 2A). A threshold was determined for each image using the 
multi-class Otsu's method in the open-source software ITK 

(Kit^\-are} [57] to separate mineralized tissue from non-mineralized 
pores. The resulting binary image was then labeled using a 3D 
connected component (CC) analysis method [41]. Objects smaller 
than 50 (im^ or larger than 1000 \lm^ were considered not to be 
lacunae and were excluded from further analysis (Fig. 3). 

Extraction of Quantitative Parameters 

The 3D image of each sample was virtually divided along the 
sample length into 3 equal-sized volumes of interest (VOI). The 
lacunae connected to the border, which could be truncated, were 

removed to avoid including bias in the analysis. Bone \-olume (BV) 
was considered as the entire mineralized tissue excluding osteocyte 
lacunae and other pores. The lacunae were segmented as 
described above and their volumes were computed for each 
VOI. The median (Lc.Vmed) a^nd the variance (LcV^ar) were 
extracted from the histogram of the lacunar volumes. Lacunar 
porosity was derived as the ratio of the total volume of all lacunae 
to the bone v()lum(; (Lc.TWBV) and th(- lacun MDD' ae density 
was defined as the number of lacunae per bone volume (N.Lc/ 
BV). Furthermore, for each VOI, the distance within which 50% 
of the mineralized bone tissue is located with respect to the closest 
lacunar surface (Lc.Distso) was computed as the median of the 
Euclidean distance transform of the bone tissue [43] . 
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Figure 3. Volume renderings of two subsections originating from the jaw (A, B) and the femur (C, D. A and C show the lacunae and 
canals from a "top-down" perspective while B and D are shown from a "front-on" perspective. The volumes of the lacunae are color- 
coded in units of |im^. 
doi:10.1371/journal.pone.0088481.g003 



Additionally, we used the reconstructed complex refractive 
index distribution, which is linearly related to the mass density, to 
compute the apparent mass density distribution (MDD"^) of each 
sample (Fig. 2B). The MDD' was calculated within the BV 
domain and was normalized by its area under the curve. Since the 
reconstructed complex refractive index might be biased due to the 
constant delta over beta ratio used in the Paganin phase retrieval 
[58], absolute values of mass density could not be retrieved and 
hereinafter the superscript r denotes relative values for all mass 
density parameters. 



However, since the different bone samples can be considered 
comparable in terms of size and composition, this allows 
quantitative comparison of the relative difference in mass density 
between the bone samples. 

Following the well-established approach for the description of 
BMDD by Roschger et al. [13], five parameters were extracted 
from the MDO"^, i.e. MDD5^|,j|jj (the mean relative mass density 
within the evaluated distribution), MDDpj,jjj^ (the most frequent 
relative mass density value), MDDlo„ and MDD^igh (the 0.5* 
and 99.5'^ percentiles), and MODp^jj^ (the full width at half 



~l 1 1 r 




mass density^ in a. u. 



Figure 4. The MDD' are shown for the three regions, with region 1, region 2, and region 3 being the top, the middle, and the lower 
sub-section of one specimen scanned, respectively. 

doi:1 0.1 371 /journal.pone.0088481 .g004 
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Table 2. Mean and standard deviation of the investigated 



properties for tine different anatomical sites and pathology are 
summarized. 






Tibia 


Femur 


Jaw 

(control) 


Jaw (BRONJ) 


Number of ROI 


9 


21 


27 


29 


Number of different 


3 


7 


9 


10 


donors 










N.Lc 


12530 


11867 


18665 


15349 


Lc.Vmed in urn' 


1 94 (27) 


224 (42) 


277 (117) 


269 (80) 


Lc.Vvar in 1 000 iim"^ 


6.4 (1.9) 


1 1 .0 (3.6) 


22.0 (11.5) 


26.6 (8.9) 


Lc.Distso in |im 


11.1 (0.4) 


15.2 (0.8) 


12.9 (1.5) 


13.9 (1.1)* 


Lc.TV/BV in % 


0.76 (0.09) 


0.45 (0.09) 


0.79 (0.32) 


0.69 (0.17) 


N.Lc/BV in 1000 mm" 


' 38(5) 


20 (2) 


27 (6) 


23 (4)* 


MDD;^„,„ina.u. 


1301 (17) 


1219 (29 


1 208 (66) 


1225 (52) 


MDDpuji^ ina.u. 


1322 (24) 


1254 (68) 


1277 (52) 


1 270 (48) 


MDDp^^,H„ ina.u. 


1 04 (38) 


149 (72) 


1 70 (66) 


142 (38) 


MDD'l„„ ina.u. 


1194 (21) 


1080 (33) 


1 048 (90) 


1073 (73) 


MDDEjig^ in a.u. 


1408 (24) 


1373 (33) 


1 360 (42) 


1 370 (45) 





Tfie * indicates that those properties are statistically significantly different 
between the BRONJ and control jaw. 
doi:l 0.1 371/journal.pone.0088481 .t002 



maxima of the distribution). The threshold of P = 0.005 was 
arbitrarily chosen and is a compromise between maintaining good 
sensitivity for low and high values in the MDD"^ and minimizing 
potential artifacts using the partial volume effect for MOD^o^,,. 

Eventually, the vascular porosity was estimated after coarsening 
the segmented volumes by a factor of five and cleaning the 
volumes from objects smaller than 1600 voxel (8575 |J,m'^). The 
canal volume (Ca.V) was quantified using voxel counting and 
canal surface (Ca.S) was determined from the number of voxel 
located within one voxel Euclidean distance to the pore boundary. 
The following parameters were quantified: ratio of canal volume 
to bone volume (Ca.V/BV), ratio of canal surface to bone volume 
(Ca.S/BV) and ratio of canal surface to canal volume (Ca.S/ 
Ca.V). 

All post-processing was done using MATLAB 2012a (The 
MathWorks Inc., Natick, MA, USA). 

Statistical Testing 

All statistical analyses were performed using the statistics 
toolbox in MATLAB. The normality of the distributions of each 
investigated parameter was determined by the Jarque-Bera test 
[59]. Differences with respect to anatomical sites, healthy and 
BRONJ groups were assessed by analyses of variance (ANOVA), 
followed by post hoc multiple comparison Bonferroni tests. The 
sample size did not allow a robust analysis of the effects of age and 
gender. All statistical results were considered significant for p< 
0.05. 

Results 

Cross sectional |iCT images and the corresponding segmented 
lacunar areas of a jaw bone sample are shown in Fig. 1. A 3D 
volume rendering of the osteocyte lacunae with subsections of one 
jaw and a femur sample image are shown in Fig. 3. These 
representative images exhibit distinct differences in the distribution 
and alignment of the lacunae in the different anatomical sites. The 



relative mass density histogram of the jaw specimen shown in Fig. 1 
is shown in Fig. 2. The normalized histograms of the three 
adjacent sub-volumes of the same specimen in Fig. 4 illustrate the 
local intra-sample variability of MDD' . 

The average properties for the three different anatomical sites 
and the BRONJ samples are summarized in Table 2 and Table 3. 

Differences between Anatomical Sites 

The average lacunar volume ranged from 194 ixm"^ in the tibia 
samples to 277 |im'' in the jaw bone. Both median values and 
variance (Figs. 5A-B) were significantiy higher in the jaw bone 
compared to the other anatomical sites. Distinct distributions of 
lacunar volume between samples from the jaw bone and those 
obtained from the peripheral skeletal sites are also illustrated in 
Fig. 6, which shows comparable, almost normal distributions for 
tibia and femur samples, but a remarkable asymmetry towards 
high volumes in the jaw bone samples. The average distance, in 
which 50% of the tissue matrix with respect to the closest lacunae 
is located (Lc.Distso), was highest in the femur, followed by the jaw 
bones and lowest in the tibiae (Fig. 5E). However, the standard 
deviation in jaw bone was considerably higher than in the other 
skeletal sites. For all investigated sites the average lacunar density 
(N.Lc/BV) was found to be larger than or equal to 20000 mm 
(Fig. 5D). Variations between the anatomical sites reflected those 
observed for Lc.Distjo, i.e. the lacunar density was highest in the 
tibia and lowest in the femur. The lacunar density in the jaw bones 
was between those of the other two sites. 

Mean, low, and high values of the relative mass density 
distribution were significantly higher in the tibia compared to jaw 
and femur (Figs. 5H^. In contrast, the heterogeneity of MDD'^ 
within the evaluated sub-volumes, as expressed by 
was higher in the jaw bone than in tibia samples. 

The average and standard deviation of Ca.V/BV, Ca.S/BV 
and Ca.S/Ca.V of all sections were found to be (7 ±5) %, 
(0.004±0.003) urn"', and (0.06±0.03) |im"', respectively. All 
values are summarized in Table 3. ANOVA revealed no 
significant differences between the different anatomical sites or 
between jaw sections from healthy donors and BRONJ. 

Differences between BRONJ and Control Jaw Bones 

Significantly lower lacunar densities (F = 5.1, p<.028) were 
observed in the BRONJ sample group (Table 2). The lower 
lacunar density was associated with higher Lc.Distso values 
(F = 6.7, p<.012). It should be noted that the values observed in 
the BRONJ samples are stiU within the range observed at other 
skeletal sites (see Table 2). All other evaluated parameters were 
statistically not significantly different when comparing BRONJ 
samples and healthy jaw bone controls. Most strikingly, the 



Table 3. Mean and standard deviation of the investigated 
morphometric properties of the segmented vessel-pores for 
the different anatomical sites and pathology are summarized. 







Tibia 


Femur 


Jaw 

(control) 


Jaw (BRONJ) 


Ca.V/BV in % 


5 (4) 


B (6) 


7 (7) 


7 (4) 


Ca.S/Ca.V in ^m" 


0.06 (0.04) 


0.06 (0.03) 


0.07 (0.04) 


0.06 (0.02) 


Ca.S/BV in ^m"' 


0.002 (0.001) 


0.003 (0.001) 


0.005 (0.004) 


0.004 (0.002) 



Ca.V: canal volume, BV: bone volume, Ca.S: canal surface. 
doi:l 0.1 371 /journal.pone.0088481 .t003 
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Figure 5. Investigated differences with respect to anatomical sites, assessed by analyses of variance (ANOVA), followed by post hoc 
multiple comparison Bonferroni tests are summarized. Significant differences between groups are indicated by a horizontal bar. If the 
significance level was reached, the p-Value and the F-Value are reported. 
doi:10.1371/journai.pone.0088481.g005 
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Figure 6. Histograms of the lacunar volumes for the three different sites are shown. Histograms are normalized to the area under the total 
number of lacunae for each site. Bin size is set to 50 )j,m^. The transparent areas indicate the standard error for each site based on the individual 
samples. 

doi:10.1371/journal.pone.0088481.g006 



lacunar volume distributions of BRONJ and control jaw bone 
samples were almost identical (Fig. 7). 

Discussion 

The present study aimed to determine the potential impact of 
alterations of the osteocyte lacunar network and mass density of 
the extracellular matrix, in the event of osteonecrosis after BP 
treatment. To address this question, we used synchrotron radiation 
phase micro-tomography with a 350 nm voxel size. We analyzed 
human cortical bone specimens extracted from the mandibular 



jaw of 10 patients suffering from BRONJ and of 9 healthy persons 
for control. In addition, we investigated samples from anatomical 
sites in which BP treatment does not usually induce necrosis, such 
as the tibia and femur. 

The imaging technique used allows the investigation of 
relatively large sample volumes in 3D, without the necessity for 
demineralization or any other tissue preparation steps (except for 
ethanol fixation and drying), and combines a large field of view, a 
very high spatial resolution and a high signal-to-noise ratio with a 
good sensitivity to mass density fluctuations [52,60]. Although the 
absolute mass density could not be derived, the derived relative 



Jaw (BRONJ) 
Jaw (control) 




1000 



Lc.V in |jm 



Figure 7. Histograms of all jaw lacunae grouped in either BRONJ or healthy bone. The shaded areas correspond to the standard error 
based on the different samples. Histograms are normalized to the absolute amount of lacunae, bin size is 50 |im^. It should be noted that even 
though the histograms of the two groups look very similar, there are differences between the histograms of individual donors. 
doi:1 0.1 371/journal.pone.0088481 .g007 
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mass density distributions (MDD^) enabled a quantitative 
comparison between the evaluated samples. 

Our results suggest that the osteocyte lacunar number, volume 
and distribution, as well as the mass density in the extracellular 
matrix, are closely linked to the anatomical site. While there are 
several studies reporting the human lacunar density based on 2D 
imaging modahties [61,62] there are only a few reporting lacunar 
density in 3D, and these are from attenuation contrast tomogra- 
phy [43] at a lower resolution (1.47 |j.m voxel size) [40,42]. The 
lacunar density values we found in the femoral samples are 
consistent with those we reported rfc(;ntly [43] and with those 
reported by Carter et al. [40] for the same anatomical region in 
women across their lifespan. In another study, significant 
variations were observed between anterior-posterior and medial- 
lateral regions with differences of up to 30% between the regions 
in single individuals, but no significant impact of age on N.Lc/BV 
[40]. Moreover, they found a significant decrease in the lacunar 
volume with respect to donor age (R^ = .46), whereas the lacunae 
from the younger group were ~30% larger (age: <50 years) than 
those of the older group (age: >50 years). In contrast, no impact of 
age on lacunar size was reported by others [27]. 

We observed pronounced differences between jaw bones in 
comparison to the two peripheral skeletal bone sites. This may be 
explained by the difiFerent origin of the cells and the difiFerent 
remodeling rates in those tissues, or the different mechanical 
environment. While osteocytes in the peripheral skeleton derive 
from the trunk lateral plate mesoderm, in the mandible they derive 
from the paraxial mesoderm [63,64]. The remodeling rate in jaw 
bone is believed to be higher than in femur and tibia [11,65], 
which leads to lower mineralization and mass density [12] in the 
jaw. The latter has been confirmed only pardy in our study by 
means of lower distribution values of MDD*^ (mean, low and high) 
and a higher intra-specimen variability (in terms of FVVHM) in the 
jaw bones compared to tibia (Fig. 5F-^. However, we observe a 
trend towards smaller values for jaw compared to femur (mean, 
low and high), which is in line with the lower remodeling rate in 
the femur reported earlier [66]. 

The majority of the osteocyte lacunar volumes were in the range 
between 100 and 400 |J.m^, which is in agreement with values 
rej)orte(l in other studies [40,42,43]. However, a small fraction of 
very large osteocyte lacunae with volumes between 400 and 
1000 nm^ was found in the jaw bone samples (Fig. 6). Most 
strikingly, this characteristic asymmetric distribution of osteocyte 
lacunar volume was almost identical in healthy jaw bone samples 
and those treated with BP (Fig. 7). The only significant effect of the 
BP treatment in jaw bone appeared to be the reduction of lacunar 
density, which was also evident in the increased average distance 
between the extracellular matrix and the closest osteocyte lacuna. 
On average, the samples treated with BP had a lacunar density 
14.8% lower than the healthy jaw bone controls, but still 15% 
higher than the femur samples. Mass density (mean, low, and high 
values) was not significantly higher in the jaws treated with BP in 
comparison to the healthy controls, although these values were 
marginally higher and the heterogeneity (FWHM) was slighdy 
lower than those in the control group. Therefore, it must be 
concluded that no remarkable changes of mass density occur in the 
extracellular matrix after BP treatment. For comparison, the mass 
density distribution values of the jaw bones treated with BP do not 
reach the values observed in tibia and femur samples. However, 
with the current analysis we cannot exclude potential local 
changes, for example directiy at the pore tissue boundary. It is well 
accepted that osteocytes actively remodel their direct extracellular 
environment [35] and it has been shown by fluorescence-labeled 
risedronate injected intravenously in a murine model that BP is 



indeed deposited around the osteocyte lacunae [67]. Interestingly, 
the results of that study have suggested that the deposition rate was 
not uniform in aU osteocytes, but decreased along with the distance 
to the next vascular channel. Therefore, a dedicated regional 
analysis of the lacunar-matrix interface merits further investiga- 
tion. Overall, except for the reduction in osteocyte lacunar 
number, the changes of the lacunar network and the mass density 
of the extracellular matrix appeared marginal. In particular, the 
BP-induced parameter alterations in the jaw stayed within the 
variations observed between different anatomical sites. Conse- 
quently, the observed parameter variations are not likely to be the 
primary causes for the development of BRONJ in the jaw. 

In this study, we have not investigated the lacunar occupation 
rate and viability of osteocytes, as the increased number of 
abandoned osteocyte lacunae is a well-documented sign of BP 
treatment in the jaw [64,68]. Moreover, the filling of empty 
lacunae with mineral has been reported in several studies [69-71], 
which is in line with the reduced lacunar density we observed in 
the jaw bones treated with BP. It is known that BPs have a high 
affinity to hydroxyapatite [63]. Therefore, we hypothesize that 
higher perfusion and turnover rates initially lead to a higher 
deposition of BP in the human jaw compared to other human sites, 
which is in line with the findings from previous studies in animals 
[11]. At higher doses BP becomes cytotoxic [72], which 
presumably promotes the gradual depletion of osteocytes and 
allows the abandoned pores to be filled with more BP-loaded 
mineral. Finally, the acidification of this tissue, caused for example 
by an inflammation, which is known to occur in BRONJ, can 
result in an excessive dissolution of BP-loaded mineral and a 
release of toxic doses of BP. Among the multifold biological factors 
promoting the development of necrosis in the jaw, the fraction of 
large lacunae (>400 |J.m'^) found in jaw bone tissue but not in the 
tissues from the other skeletal sites may be a structural factor, since 
it allows the deposition of larger amounts of highly concentrated 
BP-loaded mineral in such abandoned lacunae, which when 
washed out results in higher and thus more toxic local BP doses. 

We did not observe significant differences in vessel porosities 
between the difiFerent groups, which supports the idea that the 
observed differences in lacunar properties and mass density 
distributions brtw(-(-n the different groups are not linked to the 
potential difiFerences in vessel porosities in the investigated 
specimens. The reported pore-volume ratios are within the range 
previously shown for femoral cortical bone [43]. However, the 
quantification of vessel porosities was limited in our study, due to 
the size of the field of view, which was small in comparison to the 
average distance between individual vessels. Therefore, the field of 
view cannot be considered to be a representative volume with 
respect to the vessel network. 

This study has several hmitations. One drawback of the present 
study is that the comparison between sites and even within the jaw 
bone is hampered by our limited ability to control the exact 
anatomical location and orientation of the harvested samples and 
the underlying diagnosis. Moreover, the samples from femur, tibia, 
and jaw bones could not be collected from the same donors, and 
the duration and type of BP treatment was not uniform within the 
BRONJ group. Potential influencing factors, such as intra- 
specimen variability, age and gender, as well as the BP treatment 
conditions, could have biased our analysis. Additionally, in 
contrast to previous investigations in which we have demonstrated 
the feasibility to extract absolute mass density values from the 
phase contrast images [46,47], we report relative values in this 
study. This is due to the fact that (i) the large distance between 
sample and detector in the current configuration violated the near- 
field condition, which is an essential prerequisite for the Paganin- 
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phase retrieval used [58] and (ii) we used a constant ratio of delta 
over beta for the Paganin retrieval, even though it could have been 
different between samples with different degrees of mineralization 
[73]. 

Additionally, since we observed a bimodal distribution of the 
apparent mass densities (Fig. 4), and the samples were comparable 
small in terms of being representative for the ratio of interstitial to 
osteonal tissue, MODp^jj^ may not be an appropriate parameter to 
quantify MDD"". 

Nevertheless, the subde ultra-structural alterations observed in 
the jaws trc'at('d with BP underline the need for further 
sophisticated investigations of large tissue volume with sub-micron 
resolution, high sensitivity to local changes in mineral density and 
chemical composition of the tissue. Such changes are not likely to 
be depicted by conventional X-ray methods. Phase contrast 
tomography with voxel sizes ranging down to about 50 nm is now 
available at SR sources and may provide new hints towards the 
ultra-structural mechanisms leading to the pathogenesis of 
BRONJ. 
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